Orbital phases of fermions in an asymmetric optical ladder 
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We study a quantum ladder of interacting fermions with coupled s and p orbitals. Such a model 
describes dipolar molecules or atoms loaded into a double-well optical lattice, dipole moments being 
aligned by an external field. The two orbital components have distinct hoppings. The tunneling 
between them is equivalent to a partial Rashba spin-orbital coupling when the orbital space (s, p) is 
identified as spanned by pseudo-spin 1/2 states. A rich phase diagram, including incommensurate 
orbital density wave, pair density wave and other exotic superconducting phases, is proposed with 
bosonization analysis. In particular, superconductivity is found in the repulsive regime. 



Introduction. Orbital degree of freedom [l[ plays a 
fundamental role in understanding the unconventional 
properties in solid state materials Q. Recent experi- 
ments in optical lattices have demonstrated that orbitals 
can also be used to construct quantum emulators of ex- 
otic models beyond natural crystals. Orbital lattices are 
attracting growing interests due to their unique and fas- 
cinating properties resulting from the spatial nature of 
the degenerate states . For example, the bosonic p x + ip y 
superfluid [3H7l state has been prepared on a bipartite 
square lattice [8|, and later the other complex superfluid 
with s and p orbitals correlated was observed on a hexag- 
onal lattice @ . 

Previous study on multicomponent cold gases mainly 
focused on hyperfme states of alkali atoms HI In a 
cold gas of atoms with two approximately degenerate hy- 
perfme states, the realized pseudo-spin SU(2) symmetry 
makes it possible to emulate Fermi Hubbard model in op- 
tical lattices [l2-14|. To engineer spin-orbital couplings 
and the resulting topological phases, one has to induce 
Raman transitions between the hyp erfme states to break 
the pseudo-spin symmetry [15l-ll7|. In contrast, due to 
the spatial nature of the orbital degrees of freedom, the 
symmetry in orbital gases, such as that in p x + ip y super- 
fluid , can be controlled by simply changing the lattice 
geometry 0, 0, Gil • With a certain lattice geometry, 
a spin-orbital like coupling can naturally appear in an 
orbital gas with s and p-orbitals without Raman transi- 
tions [20|. Theoretical studies of orbital physics largely 
focusing on two or three dimensions suggest exotic orbital 
phases 
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■29( beyond the scope of spin physics. 
In this article, we study a one dimensional orbital lad- 
der with s and p orbitals coupled [2^, [3(|. We shall 
derive such an effective model for dipolar molecules or 
atoms loaded in a double-well optical lattice. The tun- 
neling rates (or effective mass) of each orbital compo- 
nent are highly tunable by changing the lattice strength. 
The coupling between s and p orbitals mimics the spin- 
orbital couplings (20[. A rich phase diagram, includ- 



ing incommensurate orbital density wave (ODW), pair 
density wave (PDW) [3ll - 433| , and other exotic supercon- 
ducting phases, is found with bosonization analysis. The 
PDW phase realized here is a superconducting phase, 
that features an oscillating Cooper pair field with a pe- 
riod of 7r. The incommensurate ODW phase has an oscil- 
lating particlc-holc pair, which tends to break the time- 
reversal symmetry. An exotic superconducting phase on 
the repulsive side is also discovered. 

Model. Consider a cold ensemble of polar molecules 



or atoms, e.g. 40 K 87 Rb [3J-l36|, 23 Na 4U K |3jj or Dy [3 



whose dipole moments are controlled by an external 
field as demonstrated in experiments. Let the ensemble 
trapped by a ladder-like optical lattice of the type stud- 
ied in [joj . That is, the lattice consists of two chains of 
potentials of unequal depth, with the molecules or atoms 
residing on the s and p orbital levels of the shallow and 
deep chains, respectively, in the tight binding regime. 
Recent experiments have reported various lattice config- 
urations of double-well potentials formed by interfering 
laser beams, showing the unprecedented tunability of the 
relative depth between two sub- wells [1, 0, [li| . The sin- 
gle particle Hamiltonian of the sp-orbital ladder is then 
given as [2(1 
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[aj(j), Q-p(j)], and a\ (a p ) is the creation 
operator for the s-orbital (p-orbital). The lattice con- 
stant is set as the length unit in this paper. The 
band structure is readily obtained by Fourier trans- 
form C 3 = J ^C(k)e ikj . The Hamiltonian in the 
momentum space reads as Hq = f ^-C' (k)H(k)C(k), 
with rl(k) = ho(k)<jQ + h{k) ■ a, where ho(k) = (t p — 
i s )cos(/c), h x (k) = 0, h y (k) = 2t sp sm(k) and h z (k) = 
— {t p + t s )cos(fc). Here Co is the identity matrix and 
o~x,y,z are Pauli matrices. The two bands are given 

by E±{k) = ho(k) ± y h y (k) + h 2 z (k), which are shown 
in FIG. [TJ The Hamiltonian is rewritten as Hq = 
I ff E p =± £: p( & ) < / > p( fc ) < / > p( fc )- Wc dcfinG an air S le li- 
able 9 by cos(#(fc)) = h z /\h\ and sin(0(fc)) = h y /\h\ to 
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save writing. Here, we only consider lower than half 
filling, i.e., less than one particle per unit cell. The 
lower band is thus partially filled and the upper band is 
empty. Since we are interested in the low-energy physics, 
the spectrum E— is linearized around the Fermi mo- 
menta ±k Fu . Here, v = A or B, and ±k FA are inner 
Fermi points and ±k FB are outer Fermi points (FIG. [I]). 
The resulting Fermi velocities are v Fv = | dE gk U=fav 
The operators capturing the low energy fluctuations 
are defined with right ('J) and left (\&) moving modes 

#a0) = <j>-(kFA+k)jJii A lk) = <f>-{-k FA +k), ttfl(fc) = 

<t>-(-k FB + k) and * B (fc) = <j)-{k FB + k). The field 
operators are introduced by ip v (x) = J ^^ v {k)e lkx and 
Tp v (x) = j *\I> u (k)e lkx . These field operators are related 
to lattice operators by 

C(j) -> \ A i> A (x)e ikFAX + \ A *4> A (x)e- tkFAX 

+ X B i) B {x)e- ikFBX + X B *$ B (x)e ikFBX , (2) 



where 



X" 



ism{0 v /2) 
cos(^/2) 



at tree level [39j. Considering strong interactions or the 
finite ranged tail of dipolar interactions, the g-ology cou- 
plings will be renormalized due to neglected irrelevant 
couplings. By manipulating the direction of dipole mo- 
ments with an external field, the interaction can be either 
repulsive or attractive [34 - 37 1 . 

We follow the notation convention of Ref. [4(| , where 
the bosonization identity takes the form 
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where r\ v is the Klein factor and d u is the dual field of 
boson field ip v . The charge and orbital boson fields are 
further introduced here by [tp c , ip ] = [(p A , (fi B ]T, with the 
matrix T given by 



T = 



V2 



1 -1 
1 1 



and their duals fields are [^cj^o] = 
Bosonized Hamiltonian density reads 



[#a,#b]T. The 



with 9 A = 8(k FA ) and B = 9{—k FB ). 

With polar molecules or atoms loaded on the sp- 
ladder, we include all momentum-independent interac- 
tions (momentum-dependent part is irrelevant in the 
Renormalization group flow |39() allowed by symmetry. 
The Hamiltonian density of the interactions is given by 

Hint = \J V ^ V ' JfJ"'] + di" JvJv' 

+53 {^aVWaV'b + ^b^a^b^a} , (3) 

where J„ =: ^tpv : and J„ =: ip*i/> v :. For the symmetric 
case t s = t p , an Umklapp process 



9u {ip* A ip A ifi B ip B + ip B il) B il)* A il) A } (4) 



becomes allowed for the reason that k FA + k FB = n. 
Since dipolar interactions between polar molecules or 
atoms decay as 1/r 3 , the leading interaction in the pro- 
posed double-well lattice setup [2fJ is 



a\(j)a s {j) 
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In the weak interacting limit, the g-ology couplings are 
related to U by gl v = U, g AB = g BA = C/sin 2 { U ^ B -), 

-,vv ff„;„ 2/a \ AB BA 



g ¥ = Usm*{e v ),g 2 A » =g 2 " A = U, 

g3 = U sm(9 A )s'm(6 B ), 

and 

g u = U cos(6U)cos(6> s ), 
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a =c/o = s]{v+ + ~grl^f -(g-rl^f , K a 
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(v FA + v FB )/2, 



and K m = J 2 J TV +9 4 , ~L , where v + 



V- = (—v FA + v FB )/2 and the transformed coupling ma- 
trices gu and g2 are given by [g] = T _1 [<?]T. The mixing 
term 'H m i x vanishes for the symmetric case with t g = t p . 

Symmetric case. For the symmetric case with t s = 
t p (FIG. [1]), the Hamiltonian has an accidental Z2 sym- 
metry, Cj — > (— \) J a x Cj and Fermi momenta are related 
by k FA = ir — k FB = k F . This Z2 symmetry implies 
that v FA = «fb, 9i AA = g BB and g 2 AA = g 2 BB ■ We 
find that the transformed coupling matrices g 2 and §4 
are diagonal and that the orbital-charge mixing term 
H m ix vanishes. In other words, the Z 2 symmetry guar- 
antees orbital-charge separation. The charge part H c 
is quadratic and the orbital part H a is a Sine-Gordon 
model (13. 

With attraction, we have K a < 1, g u > 0, and the 
Sine-Gordon term g u is relevant (flows to +00) in the 
renormalization group (RG) flow 41]. This corresponds 
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FIG. 1: Sketch of the band structure of the sp-orbital ladder. 
Red dashed lines show the level of chemical potentials, (a), 
the symmetric case with t s = t p . (b), the asymmetric case 
with t B < t p . 



to an orbital gapped phase with cos(v87r<p ) locked at 
— 1. In this phase, quantum fluctuations of ip become 
massive, and the divergent susceptibilities are the follow- 
ing: charge density wave (CDW) and PDW 31 , 32| given 
by the operators: 



-2ik FA x _ ,/,*,/, n0 2ikpB- 

oc e -^k F x e iV2^ Vc sill (V2^ ) 
Opdw(^) = ^A^ B e l{kFA+kFB}x + ^B^Ae- z{kFA+kFB)x 

Due to orbital-charge separation, the CDW and PDW 
correlation functions are readily given by 



(Ocdw(z)Ocdw(°)) K e ~" 



-2ikpx nn — K c 

j 

(Ppdw(x)4dw(0)) « (-l) x x-V K °. 



(7) 

(8) 



Since K c > 1 for attraction, the algebraic PDW order 
is dominant. In this phase, the superconducting pairing 
Osc = o-s(j)o,p(j) oscillates in space with a period of n. 

With repulsion, we have K Q > 1, and thus (73 is rel- 
evant [4l| . This gives an orbital gapped phase with 
cos(-\/87r# ) locked at 1, because 33 < 0. The fluctua- 
tions of i? are massive, and the divergent susceptibilities 
are ODW and superconducting SC + given by the opera- 
tors: 



°SC+{ X ) = 1pAll>A + 1pBi>B 

oc e~ iV ^^ cos(v&„) 

Since K c < 1 for repulsion, the dominant algebraic order 
here is ODW, for which the correlation function is given 
by 



(OoD\v(a;)Oj )DW (0)) oc e 



-i(k F A—k F B)x„—K c 



(9) 



In the ODW phase, the particle-hole pairing in terms 
of lattice operators reads OqdwC?) = Cj°yCj 



ODW order is incommensurate with an oscillation pe- 
riod 2-Kj(kpA — kps) in real space. If we go beyond the 
one-dimensional limit and consider small transverse tun- 
nelings [2(J, a true long-range ODW order (Oodw^)) oc 
e i(k FA -k FB )x j g ex p ec ted. Such an order breaks time- 
reversal symmetry. 

The ODW and PDW phases predicted by Bosoniza- 
tion analysis are further verified in numerical simula- 
tions with matrix products states. The superconducting 
correlation C S c(f - j) = {ap{j)a\{j)a s {j')a p {j')) and 
the orbital density wave correlation Codw(j' — j) = 
(Cj(jyCjCj,<TyCji) are calculated. FIG. [5] shows the 
Fourier transform of these correlations, defined by C(k) = 
YljM @(j) e ~ ■ The sharp peaks of Csc{k) at momenta 
±7r on the attractive side tell the quantum state has a 
PDW order shown in Eq. [8] On the repulsive side sharp 
dips of Codw(^') at finite momenta verify the incommen- 
surate ODW order shown in Eq. [9] With numerical cal- 
culations, we also find the existence of PDW phase in the 
strongly attractive regime if t s ^ t p . 
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FIG. 2: The phase diagram of the symmetric sp-orbital ladder 
with t s = t p = 2t sp . Codw(A:) and Csc(fc) show the Fourier 
transform of the orbital density wave and superconducting 
correlations, respectively. The numerical results are calcu- 
lated for a system of 70 particles in a ladder of size 100. In 
the upper (lower) inset, the interaction U = 3t 3 (U — — 3t 3 ). 



Asymmetric case. For the asymmetric case — t s < t p 
(FIG.fl]), the Fermi velocity vps > vpa and the orbital- 
charge separation no longer holds. Thus, the orbital and 
charge degrees of freedom cannot be treated separately. 
The other difference with the symmetric case is that the 
Umklapp process g u does not exist. Since the effects of 04 
couplings are just to renormalize the Fermi velocities [42- 
[45| | . For simplicity, we do not consider such effects and 
set 54" = here. The one-loop RG equations are given 
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by E3, 



dg 



V V 



dl 

dm 
dl 
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VFv 



(10) 



where I is the flow parameter (/ — > oo) and P = A (£?) 
for v = B (A). The RG flow of the Sine-Gordon term g 3 
is obtained as 



\C\/g 3 (l) 



F 



-^( 93 Y)J 2 -^l + F-i 



\O\/g 3 (0)} (11 
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VFAV FB+v'i 
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(12) 



(13) 



(14) 



The function i* 1 is the hyperbolic function "sinh" (the 
trigonometric function "sin") if C > (C < 0). When 
y > 0, 33 always flows to 00 and the system is in some 
gapped phase. When Y < 0, g 3 flows to 00 only if C < 0. 
<? 3 is irrelevant only if C > and y < 0. In the weak 
interacting regime, Y/U = ^~ 



2v FA 



2v FB 



We 



will consider the regime Y/U > (this condition holds 
when t sp is weak compared with t s + t p ) in the following. 

With repulsion (U > 0, Y > 0), 33 is relevant and 
flows to —00 in RG flow. Then the dual orbital field 
$0 is locked with cos(v / 87n9 ) = 1 and its fluctuations 
$ are massive. The key effect of orbital-charge mixing 
can be seen from its modification of the dynamics of the 
conjugate fields, given as 

K„ „ . K a u Tl 



n« 



n 



-dtdo 



K, 



1 



-d t (p c 



d x tp 



dx'&o, 



(15) 
(16) 



u c K c K c u c 

where II ro is the conjugate field of vj. The Lagrangian is 
constructed by £(x,t) = He o d t o + U iPc d t ip c ~ H- With 
massive fluctuations of i) integrated out, the Lagrangian 
of the charge field tp c is given by 



1 

2^ 



(d t ip c ) 2 - u(d x (p c y 



O ((d<p c ) 



(17) 



with the renormalized Luttingcr parameter and sound 
velocity given by 



7 



K C K„ 



u c K c K a 



(18) 



(19) 



To zeroth order in the interaction U, the renormalized 

-1/2 



Luttingcr parameter is 7 
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Our result 



reproduces the pcrturbative result [45( when the orbital- 
charge mixing term is small. The diverging susceptibili- 
ties are ODW and SC + , and the corresponding correla- 
tion functions are given as 

(O sc+ (x)c4 + (0)) ex x- 1 '^ (20) 
(Oodw(z)OUv(O)) « e -^ A - k -^x-\ (21) 

With sufficiently weak repulsion 7 > 1, the dominant 
order is SC + , of which the pairing in terms of lattice op- 
erators is Osc = o, s (j)a p (j). We emphasize here that 
this pairing does not oscillate in real space. Such a su- 
perconducting phase arises in the repulsive regime due 
to the orbital-charge mixing and the pinning effect of the 
dual orbital field # . The Sine-Gordon term 173 caus- 
ing this pinning effect is finite only when the coupling 
of sp-orbitals t sp is finite, and g% is monotonically in- 
creasing when t sp is increased. Thus the transition tem- 
perature of this repulsive superconducting phase can be 
increased by tuning t sp , which makes this exotic super- 
conducting phase potentially realizable in experiments. 
With stronger repulsion, the renormalized Luttinger pa- 
rameter 7 decreases. Eventually with repulsion larger 
than some critical strength, we have 7 < 1, and the re- 
pulsive superconducting phase gives way to the ODW 
phase. Our numerical results cannot distinguish the 
ODW and SC phases in the weak interacting regime due 
to computational complexity j4(| . 

With attractive interaction, the condition Y/U > 
gives y < 0. Thus g 3 is relevant and flows to +00 when 
C < 0. The Sine-Gordon term cos(v / 8~7ri9 ) is locked 
at —1, and the dominant order is superconducting SC _ , 
given by 



O, 



sc- 



IpAlpA ~ 



sin(\/27n? ). 



(22) 



In numerical simulations wc find the SC _ phase compet- 
ing with PDW in the strongly attractive regime. When 
cfy is irrelevant (C > 0, Y < 0), the orbital ladder is in 
a two component Luttingcr liquid phase exhibiting two 
gaplcss normal modes and each mode is a mixture of or- 
bital and charge. 
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